Modelling non-dust fluids in cosmology 
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Often in cosmology Newtonian physics is used to model the dynamics of matter inhomogeneities. 
For pressureless dark matter, or 'dust', this approach gives the correct results, but for scenarios 
in which the fluid has pressure this is no longer the case. In this article, we explicitly highlight 
the relationship between the variables in Newtonian and cosmological perturbation theory, showing 
exact equivalence for pressureless matter, and giving the relativistic corrections for matter with 
pressure. As an example, we focus on the scalar field dark matter model, which has recently gained 
popularity but which has non-zero pressure perturbations. We discuss some problems which may 
arise when modelling this theory with numerical simulations, and when using CMB Boltzmann 
codes. 

PACS numbers: 98.80.Cq |arXiv:1207."l870| 



I. INTRODUCTION 

Current observations indicate that the universe in 
which we live is almost homogeneous and isotropic. How- 
ever, it is also known that small initial departures from 
homogeneity and isotropy give rise to the structures we 
observe today. Thus, while the universe is well approxi- 
mated on large scales by a homogeneous and isotropic 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space- 
time, the existence of large scale structure and inhomo- 
geneities in the Cosmic Microwave Background (CMB) 
tell us that this is just an approximation pQ. 

In order to model inhomogeneities we utilise a tech- 
nique that is well established in many branches of physics 
and applied mathematics, namely, we take an approxi- 
mate solution and add small perturbations. In cosmol- 
ogy, this process is called cosmological perturbation the- 
ory, and requires the addition of small inhomogeneous 
perturbations on top of the FLRW background, such 
that the system still solves Einstein's field equations (e.g. 
Refs. [2H1|)- However, when considering dynamics on suf- 
ficiently small scales, and for fluids which are pressure- 
less, it is sufficient to use Newtonian physics [5;. Inho- 
mogeneous perturbations of Newtonian cosmology have 
been studied for many years and these are the equations 
that are used when performing large numerical simula- 
tions of galaxy formation [6! . However, as simulations get 
more sophisticated and the observations more accurate, it 
is important to consider the contribution of those compo- 
nents which exhibit pressure perturbations. For example, 
to fully account for the effects of inhomogeneous scalar 



fields as dark energy (e.g. [7]) or dark matter [5HTT] in the 
formation of structures, one should employ more general 
equations with the input of general relativity. 

In this paper we re-visit the question of how to re- 
late Newtonian and cosmological perturbation theories. 
We show that for pressureless systems the equations gov- 
erning cosmological perturbations of an FLRW space- 
time reduce to the equivalent Newtonian equations on us- 
ing gauge invariant variables defined in different gauges, 
namely the metric potential in the longitudinal gauge, 
and the density contrast in the comoving gauge |12) . 
Drawing on this equivalence, we then investigate the situ- 
ation for fluids with pressure. We find that one can write 
the Poisson equation in the usual way, but that the con- 
tinuity and Euler equations then differ, depending upon 
the equation of state parameter and the adiabatic sound 
speed. 

We then go on to study how to relate the two per- 
turbation theories in a scalar field model. We find that, 
again, the Poisson equation is identical to the Newto- 
nian case, but that the Euler and continuity equations 
differ, depending now on the equation of state parameter 
and the effective speed of propagation of perturbations 
through the system. Finally we discuss the Jeans scale 
in the scalar field dark matter models, where the back- 
ground equation of state parameter is zero. We find that 
the scales depend upon Cp h , the phase speed, or speed of 
propagation of perturbations. We close, in Section |IV| 
with a brief discussion. 



II. MODELLING INHOMOGENEITIES 



A. Newtonian perturbations 
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Let us first study the theory of perturbations in Newto- 
nian physics. We consider inhomogeneous perturbations 
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about a homogeneous background, and so the energy den- 
sity p is 

p(x,t) = p(t)(l + S(x,t)), (2.1) 

where p(t) is the homogeneous background energy den- 
sity and S(x,t) is the inhomogeneous density contrast. 
On introducing the inhomogeneous Newtonian potential, 
$n(^, t) and fluid velocity v(x,t), the linearised conser- 
vation and Euler equations are, respectively, ITS] 



1, 



-V- v = 0, 
a 

il+ Hv = --V$ N - -4v£P. 
a ap 



(2.2) 
(2.3) 



where H — a/ a is the Hubble parameter, a dot denotes 
a derivative with respect to coordinate time and SP de- 
notes the pressure perturbation. The Newtonian poten- 
tial and the density contrast are then related through the 
Poisson equation 



V 2 $ N = 4:nGa 2 p6 . 



(2.4) 



where the Laplacian is defined as V 2 = V • V. On utilis- 
ing the relationship between the energy density perturba- 
tion and the pressure perturbation for a barotropic fluid, 
SP = c 2 5p, we can combine the fluid equations into a 
second order equation 



(PS 
dt 2 



8S 

2H— = 4irGpS + c 2 V 2 S . 



(2.5) 



B. Cosmological perturbations 



While the theory of Newtonian perturbations is suf- 
ficient for modelling small scale physics involving only 
pressureless dark matter particles, the dynamics of the 
universe are governed by general relativity. Therefore, 
we must consider relativistic perturbation theory. Since 
Einstein's theory relates the geometry of the universe to 
its matter content, we must consider perturbations of 
both the matter and the FLRW spacetime metric. 

The most general, linear scalar perturbations to the 
FLRW metric are [1 H HI 

ds 2 = a 2 (r/) - (1 + 2<f))dr] 2 + 25 l dx l d'q 

+ {(1 - 2i>)5 ij + E^}dx l dxA , (2.6) 

where we now use conformal time ry, related to coordinate 
time t through dt = adrj. A unique problem which arises 
in the relativistic theory is the problem of gauge depen- 
dence. Since general relativity is covariant, and splitting 
the spacetime into a background and a perturbation is 
not a covariant process, we introduce extra spurious co- 
ordinate dependence (see, e.g., Refs. [HE3HE])- This can 
be resolved in a systematic manner, as was first shown 



by Bardeen |12) , by considering gauge-invariant variables 
- quantities that do not change under a gauge transfor- 
mation. A popular choice of variables amount to setting 
E and B to zero, resulting in the FLRW metric in the so- 
called longitudinal or Newtonian gauge, with the gauge 
invariant variables $ and ^ 

ds 2 = a 2 {i 1 )\-{l + 2^)dr j 2 + {l-2^)5 ij dx l dxA . (2.7) 

The conservation and Euler equations are then reduced 
from the form without fixing the gauge 



Spf + m{5p + SP) - 3( j o + P)tp' 

+ (p + P)W 2 (v + E r ) =0 



V 



'+(1-3^)™ 



(SP+ 2 V 2 IL 
p + P \ 3 



(2.8) 

= 
(2.9) 



to the following expressions, where we neglect anisotropic 
stresses, II, 

^ + (p + P)(V 2 i; £ -3*') = -3H(5 Pe + 5P e ), (2.10) 

[{p + P)v t ] ' + m{p + P)v t + (p + P)$ + SP i = 0. 

(2.11) 

Here the subscript I denotes matter variables in the 
longitudinal gauge and v is the velocity potential, i.e. 
v l = V l v with v % the fluid three velocity, and V = v + B. 

We can then specialise to a barotropic fluid with equa- 
tion of state P — w(r))p, and whose pressure perturba- 
tion can be related to the energy density perturbation 
through 



SP = c 2 5p, 



(2.12) 



where c 2 = P'/p' is the adiabatic sound speed 1 . Thus, 
for this system, Eqs. (2.101 and (2.11) become 



S' t + (1 + w)(V 2 v e - 3*') = m{w - c 2 )S e , (2.14) 



v , i +n{l~3c 2 s )v e + ^ 



1 + w 



-Si = 0. 



(2.15) 



The Einstein equations then give that $ = $ (in the 
case of zero anisotropic stress, as is true for any perfect 
fluid), and the Poisson equation is 



V 2 * = AiTGa 2 p 



Si — 3H(l + w)vi 



(2.16) 



1 The adiabatic sound speed, equation of state parameter and its 
derivative are related through 



1 + w 



-3H(c 2 s - w) . 



(2.13) 



For the case where w is constant, these relationships guarantee 
that w = c?. 
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In order to draw an equivalence between this and the 
Newtonian Poisson equation, we can set ^ = <&n- Fur- 
thermore, on considering that the energy density per- 
turbation transforms under the gauge transformation 



Srj as 



Sp = Sp + p'Srj , 



(2.17) 



we get for the comoving density contrast, in terms of the 
longitudinal density contrast and the velocity perturba- 
tion P3], 

5 c = 5 l -m{l+w)vt. (2.18) 
Then wc obtain the Poisson equation [T7] 



V 2 $ N = Ai:Ga 2 p5 c 



(2.19) 



which is equivalent to Eq. (2.4) 



C. Correspondences 

In this Section we relate the the relativistic equations 
to the Newtonian equations. We take heed from the Pois- 
son equation which, as stated above, relates the density 
contrast in the comoving gauge to the metric potential 
in the longitudinal gauge: 



V 2 * = ^Ga 2 p8 c 



(2.20) 



In order to obtain an equivalence for the set of equa- 
tions, we must first ensure that we are consistent with 
the density contrast that we use. With the aid of the 
background Friedmann equations and one of Einstein's 
perturbed equations, 



= -HV - 47rGa 2 (l + w)pv e . 



(2.21) 



we can rewrite Eq. (2.10) as 



S' e + m{c 2 s - iv)5 e + (1 + to) (SUV + V 2 v e 



0. 



(2.22) 



Using now the relationship between Si and S c , Eq. ( |2.18 1, 
and its derivatives, the continuity equation can be writ- 
ten as 



5c - 3Hw5 c + (1 + w)W 2 v e = . 



Furthermore, the Euler equation, (2.11), is 
v', +Hv e = - 



1 



S r - 



(2.23) 



(2.24) 



From this we can see that, in the case of a pressure- 



less dust for which w 
reduce to 



0, the evolution equations 



V V = , 
■Uvi + * = , 



(2.25) 
(2.26) 



which can be recast as a second order differential equa- 
tion as 



S" + HS C = AiiGa 2 pS c 



(2.27) 



which governs the evolution of matter density perturba- 
tions. Here the second term is a suppression of the per- 
turbations with the expansion of the universe, and the 
term on the right hand side sources the growth of pertur- 
bations due to the gravitational instability. This is the 
equivalent form to the Newtonian perturbation theory 
(as expected for a non-relativistic species), as obtained 
in Eq. (2.5). This equivalence was discussed in greater 



detail in Ref. [15] . 

However, should one wish to consider fluids other than 
dust, with non-zero pressure perturbations, then the rel- 
ativistic equations must be used. This is especially im- 
portant for, say, hot dark matter, or for a system con- 
taining dark energy perturbations. In the following, we 
will consider the case of scalar field dark matter models. 



III. SCALAR FIELD DARK MATTER 

A model of dark matter which is receiving increased 
attention is the scalar field dark matter model (SFDM), 
specifically in the form of a Bose-Einstein condensate 
[H?l |2"U] . the theory is robust enough to study the nu- 
merical problem of structure formation (recent attempts 
are [2J-23J). The basic model features a canonical scalar 
field field oscillating at the bottom of its potential with a 
period much smaller than the Hubble time and any other 
dynamical times. In consequence the SFDM can be re- 
garded as a fluid with zero effective pressure. However, 
in contrast to cold dark matter (CDM), this component 
exhibits pressure perturbations. 

It is well known that a scalar field system cannot be 
modelled as a barotropic fluid (except on super-horizon 
scales ;24J), and in fact a consistent fluid equivalence 
must consider the more general perfect fluid |25j . In- 
deed, as elucidated in Ref. [35], in order to interpret a 
scalar field as a fluid, the equation of state parameter, 
w and adiabatic sound speed c 2 differ. Furthermore, the 
speed with which pressure perturbations propagate is de- 
scribed by the effective sound speed, or phase speed, c 2 h 
defined, in the fluid rest frame, as 



SP 
Jp 



(3.1) 



For a scalar field with a canonical kinetic term, this is 
equal to unity. However, for a non-canonical scalar field 
with pressure and energy density depending upon both 
the field, ip, and its kinetic term, X = —^g^d^ipd^ip, 
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this can differ from one and is given by [2"41, 127] 2 

P.x 



r 2 _ 



(3.2) 



For SFDM, the pressure perturbation is no longer pro- 
portional to the energy density perturbation. Instead the 
relationship between pressure perturbations and energy 
density perturbations is 



SP = c 2 Sp + (c 2 - c 2 s ) Sp + p'(v + B) 



(3.3) 



The second term on the right hand side is often referred 
to as the non-adiabatic pressure perturbation 5P na &. In 
this case, Eqs. (2.10) and (2.11) become 



<5; + (l + u;)(V 2 W£-3*') 

= m{w - c? vh )8 t + m 2 {l + w)(c 2 h - c 2 )v e , (3.4) 



v' t + H(l - 3w)v e + $ + —^-5i + - 

1 + w 1 

C ph 



-V£ 



3H(c 2 h - c 2 )ve . 



(3.5) 



As shown previously, in order to write the equations 
in the variables employed in simulations, we introduce 
the density contrast in the comoving gauge. Following a 
similar procedure, using the background and perturbed 
Einstein equations, we obtain 

S' e + 3H(c 2 ph - w)Si + 3H 2 (1 + + (1 + w)V 2 Vi 

1 



+ m z (i + w)(ci - c ; h + 



-)v t = 0, 



P h ' 2 
which can then be written as 

S' c - 3Hw6 c + (1 + w)V 2 v e = . 



(3.6) 



(3.7) 



Interestingly, this is identical to the equation for the 
barotropic fluid presented earlier. However, the Eulcr 
equation, Eq. (2.24), does not take the same form, and 



can be written in terms of the comoving density contrast 
as 



1 



Sr. 



W 



(3.8) 



The combination of these two equations yields a Klein- 



Gordon equation - a generalisation of Eq. (2.27) for S 
with an extra Laplacian term, 



cl h V 2 S c +H(3w-l)S' c 



m. 1 



3w + -w ) S c 



0. (3.9) 



From Eq. (3.9) one can determine the instability scale 



for density perturbations. For the most general, scalar 
field case we obtain the Jeans wavenumber 



3 n 1 



3w — —w 



(3.10) 



This scale has been associated to a power spectrum cut- 
off, and is a characteristic feature of non-dust compo- 
nents. In the case of the scalar field acting as a dark 
matter component, for which w = 0, the Jeans' wave- 
length reduces to 



A, 




(3.11) 



This is equivalent to the field calculation (see, e.g., 
Ref. [35]) but with the subtle difference that the effec- 
tive sound speed, c p h, takes the place of the adiabatic 
sound speed c s (see also [1H])- This corresponds to the 
fact that the speed of propagation of perturbations in a 
scalar field system is given by c 2 h . The effective sound 
speed for each canonical SFDM model is equivalent and 
equal to 1, though for non-canonical models this differs 
(see e.g. [3D]). The specific scales for each model, as well 
as the growth of perturbations will be treated elsewhere. 

These simple results show the importance of consider- 
ing the system with no approximations and argue for the 
use of the system of equations ( 3.7 1 and ( 3.8 ) in the forth- 



coming simulations of structure formation in the SFDM 
model. 

Another important comment is that care must be 
taken when studying the SFDM model using CMB codes. 
Many of the popular Boltzmann codes are written in the 
synchronous gauge, a gauge specified by demanding that 
cj> = = B. However, as has been known for some time, 
these conditions alone do not fix the gauge. That is, 
the synchronous conditions alone are not a proper gauge 
choice, and one needs an additional condition. The con- 
dition usually chosen, and which is employed by, e.g., 
CAMB [3T], is to set the velocity perturbation of the 
cold dark matter to zero. 

This can be done since the dark matter is assumed to 
have zero pressure perturbations [32) . But this is not 
possible for models with pressure perturbations, since 
imposing this condition will result in inconsistencies in 
the theory. Consequently, in order to use CMB codes to 
study the SFDM model, one must either include an ex- 
tra CDM component to allow for this condition - which 
arguably reduces the value of the theory - or use a code 
that is not written in the synchronous gauge. 



IV. DISCUSSION 



2 Note that the equivalence between Eqs. ( |3.l| and j3.2| is only 
true for a single scalar field, which is the case we study in this 
paper. 



In this paper we have revisited the issue of relating 
Newtonian perturbation theory to relativistic perturba- 
tion theory in cosmology. Having reviewed both pertur- 
bation theories, we then explicitly showed how one relates 
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the two for the case of a pressureless fluid, such as dark 
matter. As is quite well documented in the literature on 
cosmological perturbation theory, in order to relate the 
two approaches one must use the Newtonian gauge with 
the comoving energy density perturbation. 

We then considered the extension of this for fluids with 
pressure, taking our lead from the dark matter case, and 
using the Newtonian gauge with the comoving energy 
density perturbation, which results in the Poisson equa- 
tion being identical to the Newtonian case. In this case 
the continuity equation differs, depending upon the equa- 
tion of state parameter, and the Euler equation depends 
upon Cg. 

The major application that we explored in Section [TTT| 
regards scalar field dark matter. In this case, the dark 
matter species is no longer pressureless, but is instead a 
scalar field with canonical kinetic term. In addition to the 
adiabatic sound speed and equation of state parameter, 
which are introduced for a general perfect fluid, a scalar 
field also has a speed of propagation of perturbations, 
which we dubbed the phase, or effective sound, speed. 
The evolution equations for this scalar field fluid in the 
longitudinal gauge depend on all the parameters. How- 
ever, on writing the equations in terms of the comoving 
energy density perturbation, the Poisson and continuity 
equations reduce to those of perfect fluid form. The Eu- 
ler equation, on the other hand, still depends upon all 
the parameters. Thus we have shown that, when study- 
ing the scalar field dark matter model and treating it 



as a fluid, one cannot simply use the pressureless cold 
dark matter equations without potentially finding erro- 
neous results. Instead, one must use the equations given 
m Section Hill 

This non-equivalence of the relativistic theory to the 
Newtonian theory is also important in systems containing 
more than one fluid. For example, as shown in Ref. [33], 
in a system containing normal CDM and a dark energy 
component, Eq. (2.271 is no longer satisfied, since one 



must take into account the dark energy perturbations. 
This may also be important when attempting to model 
a system with dark matter and another species. 
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